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A  mathematical  model  is  formulated  to  describe  incipient  breaching  of  coastal  barrier 
islands.  The  model  is  based  on  the  assumptions  of  idealized  breach  morphology  and  is  in¬ 
tended  to  describe  the  growth  of  breaches  prior  to  possible  closure  by  longshore  sediment 
transport.  The  two  coupled,  non-linear  equations  governing  breach  width  and  depth  are 
solved  analytically  for  special  cases.  The  analytical  solutions  explicitly  exhibit  an  expo¬ 
nential  behavior  in  breach  dimensions  and  reveal  that  the  macroscale  process  of  breach 
growth  is  controlled  by  seven  variables:  initial  width  and  depth  of  the  breach,  equilibrium 
width  and  depth  of  the  breach,  width  of  the  barrier  island,  and  maximum  or  initial  net 
sediment  transport  rates  at  the  bottom  and  sides  of  the  breach.  The  literature  is  reviewed 
to  compile  general  properties  of  coastal  breaches,  and  sensitivity  testing  shows  the  model 
to  be  compatible  with  those  observations.  The  model  is  applied  to  simulate  the  1980  breach 
at  Moriches  Inlet,  New  York,  and  reasonable  agreement  is  found. 

Keywords’.  Breach;  inlet;  barrier  island;  mathematical  model;  Moriches  Inlet;  morphology 
modeling. 


1.  Introduction 

A  coastal  breach  is  a  new  opening  in  a  narrow  landmass  such  as  a  barrier  island 
or  barrier  spit  that  allows  water  to  flow  between  the  water  bodies  on  each  side. 
Breaching  is  a  common  natural  occurrence  on  barrier  islands  and  spits,  and  at 
ephemeral  river  mouths  and  coastal  lagoons  and  ponds.  Breaches  are  also  induced 
artificially  for  flood  protection,  to  increase  water  quality,  and  as  an  environmental 
enhancement.  Natural  breaching  can  be  initiated  in  three  ways:  (1)  inundation  from 
the  sea  combined  with  wave  action,  (2)  elevated  water  level  in  the  back  bay  or  river, 
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causing  piping  and  liquifaction,  and  (3)  narrowing  of  a  barrier  island  or  spit  by 
reduction  of  sediment  supplied  through  longshore  transport  (Pierce,  1970;  Kraus, 
Militello,  and  Todoroff,  2002).  Coastal  barrier  breaching  is  expected  to  become  more 
prevalent  with  rise  in  sea  level,  erosion  of  the  coast,  and  continued  length  of  service 
of  jetties. 

Once  breaching  occurs,  tidal  exchange  and  river  discharges  will  typically  widen 
the  initial  opening.  A  breach  may  close  naturally,  but,  if  the  tidal  exchange  is  strong 
and  longshore  sediment  transport  weak,  it  can  increase  in  size  and  become  a  new 
inlet.  An  inlet  created  by  breaching  may  compete  for  stability  with  existing  inlets 
in  the  same  bay  system,  promoting  their  closure.  A  breach  adjacent  to  a  jetty  has 
potential  for  undermining  and  isolating  the  structure,  and  it  will  convey  a  portion 
of  the  flow  that  would  otherwise  scour  the  navigation  channel,  promoting  shoaling 
of  the  channel.  In  addition,  nearshore  and  beach  sediment  that  had  been  protecting 
the  structure  and  the  shore  may  move  through  the  breach  and  into  the  navigation 
channel,  increasing  dredging  requirements  and  adding  unanticipated  cost  to  inlet 
entrance  maintenance.  A  trend  toward  closure  of  an  existing  inlet  will  make  the 
navigation  channel  unreliable,  as  well  as  alter  the  environment  because  of  changes 
in  water  level,  circulation,  and  salinity. 

Breaches  can  cause  loss  of  human  life,  property,  infrastructure,  and  transporta¬ 
tion  corridors;  endanger  navigation  and  stability  of  adjacent  inlets;  and  harm  the 
environment  by  the  exposure  of  bay  perimeter  to  sea  waves  and  by  the  change  in 
bay  salinity.  The  cost  of  breach  closure  can  be  high,  so  unintended  breaches  that 
will  be  a  concern  should  be  avoided  by  counter  measures  and  breach  contingency 
planning.  Most  breaches  open  rapidly  and  then  gradually  evolve  over  a  period  of 
weeks  to  months  if  they  continue  to  remain  open.  Quantitative  predictive  tools  are 
necessary  to  assess  vulnerability  of  coastal  barriers,  design  breach-prevention  mea¬ 
sures,  develop  breach-closure  plans,  estimate  the  fate  of  a  breach,  and  evaluate  the 
consequences  of  a  breach  to  the  neighboring  inlets,  beach,  and  estuarine  system. 
Measurements  and  reliable  predictive  models  of  coastal  breaching  are  lacking. 

This  paper  introduces  a  mathematical  model  of  the  breaching  of  alluvial  coastal 
barriers  based  on  a  macro-scale  or  morphologic  approach.  The  model  provides  a 
heuristic  framework  for  understanding  the  breaching  process  and  makes  apparent 
key  dependencies  that  control  the  growth  of  a  coastal  breach. 


2.  Breaching  Process 

Review  of  the  literature  indicates  that  little  quantitative  information  is  available 
on  the  process  of  coastal  barrier  breaching  and  that  predictive  capability  is  lacking 
(Basco  and  Shin,  1999;  Kraus,  Militello,  and  Todoroff,  2002;  Kraus  and  Wamsley, 
2003).  During  hurricanes  and  storms,  large  numbers  of  ephemeral  breaches  can  open 
(e.g.  Texas  coast,  USA  —  Price,  1963;  Hayes,  1967;  and  Pierce,  1970;  Louisiana 
coast,  USA  —  Wright,  Swaye,  and  Coleman,  1970;  and  New  Brunswick,  Canada 
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—  Greenwood  and  Keay,  1979).  In  a  storm  report  documenting  the  damage  of  the 
September  1967  Hurricane  Beulah,  the  US  Army  Engineer  District,  Galveston  (1968) 
states  that  an  inspection  team  found  31  breaches  along  48  km  of  North  Padre  Island, 
Texas,  barrier  beach.  A  book  edited  by  Leatherman  (1981)  collects  classic,  primarily 
descriptive  papers  on  the  processes,  sedimentation,  and  morphology  of  overwash  and 
breaching.  Other  papers  are  discussed  below  to  supplement  and  update  information 
in  Leatherman  (1981). 

Dent  (1935)  documents  the  breaching  of  a  90-km  long  barrier  beach  at  Ocean 
City,  Maryland,  during  an  August  1933  hurricane.  This  breach  was  quickly  stabi¬ 
lized  with  jetties  to  form  Ocean  City  Inlet  and  protect  its  federal  navigation  entrance 
channel.  Terich  and  Komar  (1974)  document  the  breaching  of  Tillamook  Spit,  Ore¬ 
gon,  that  experienced  a  reduction  in  longshore  sediment  source  by  the  construction 
of  a  jetty.  At  high  tide,  the  water-filled  breach  at  Tillamook  was  nearly  1.2  km 
wide.  It  was  closed  four  years  later  by  the  construction  of  a  setback  rubble  stone 
revetment,  with  a  sandy  beach  filling  the  indentation. 

Rice  (1974)  describes  closure  conditions  for  the  mouth  of  the  Russian  River,  Cal¬ 
ifornia.  He  noted  an  apparent  causal  relation  between  opening  of  the  river  mouth 
during  strong  stream  flow,  and  closure  of  the  mouth  to  the  season  of  largest  waves 
and  greatest  longshore  transport.  Nishimura  and  Lau  (1979)  discuss  possible  means 
of  minimizing  closure  of  the  mouths  of  small  streams  and  propose  innovative  struc¬ 
tures  to  promote  hydraulic  breaching  by  early-arriving  storm  flows  as  the  eroding 
agent.  Kraus,  Militello,  and  Todoroff  (2002)  document  breaching  of  a  barrier  spit 
by  elevated  water  level  in  Stone  Lagoon,  California,  and  discuss  the  common  occur¬ 
rence  of  breaching  of  the  Humboldt  State  Park,  California,  lagoons  at  the  end  of 
spring,  the  season  of  maximum  precipitation.  The  barrier  spits  fully  enclosing  the 
Humboldt  lagoons  breach  from  the  lagoon  side  when  the  water  level  rises  to  about  4 
m  above  mean  sea  level.  The  lagoon  breaches  open  to  maximum  dimensions  within 
several  hours,  similar  to  dike  failure,  remain  open  at  low  tide  for  days  to  weeks  but 
may  be  passable  at  high  tide,  then  close  through  infilling  by  longshore  sediment 
transport. 

Smith  and  Zarillo  (1988)  document  the  growth  of  an  artificially  induced  breach  at 
Mecox  Pond,  Long  Island,  New  York,  that  remained  open  during  September  10-18, 
1985.  Mecox  Pond  and  neighboring  coastal  ponds  are  breached  by  digging  a  channel 
to  the  Atlantic  Ocean  when  precipitation  raises  the  level  of  the  ponds  to  threaten 
flooding  of  neighboring  property.  Figure  1  shows  (a)  the  west  bank  of  the  channel, 
and  (b)  the  length  of  a  breach  at  Mecox  Pond  that  was  cut  from  a  narrow  pilot 
channel  in  February  1998.  Steepness  of  the  walls  of  the  highly  rectilinear  breach  is 
evident. 

Terchunian  and  Merkert  (1995)  document  two  barrier  island  breaches  formed 
at  Westhampton,  Long  Island,  New  York,  for  which  the  first  and  initially  larger 
breach  (Fig.  2)  closed  by  longshore  transport  and  moderate  engineering  actions.  A 


514  N.  C.  Kraus 


(a)  Looking  west  at  breach  cut  through  beach  berm  at  Mecox  Pond;  Atlantic  Ocean  to  left. 


(b)  Looking  north  from  Atlantic  Ocean  to  Mecox  Pond. 

Fig.  1.  Mecox  Pond,  breached  on  morning  of  February  14,  1998;  photographs  taken  in  afternoon  of 
February  14. 
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Fig.  2.  Pikes  Inlet,  New  York,  December  18,  1992,  about  6  days  after  inception  of  breaching. 


neighboring,  initially  smaller  breach  widened  because  of  updrift  trapping  of  sediment 
by  a  groin  field.  The  latter  inlet,  called  Little  Pikes  Inlet,  grew  from  approximately 
30  m  to  1.5  km,  more  than  five  times  the  width  of  the  nearby  Moriches  inlet  to  the 
west  that  is  stabilized  by  dual  jetties. 

General  observations  and  a  few  case  studies  of  breaching  that  occurred  next  to 
jetties  have  also  been  made,  such  as  for  the  January  1980  breach  at  Moriches  Inlet, 
New  York  (Schmeltz  et  al. ,  1982;  Sorensen  and  Smeltz,  1982)  and  the  December 
1993  breach  at  Grays  Harbor,  Washington  (Arden,  2003).  The  breach  at  Grays 
Harbor  (Fig.  3)  was  strongly  flood  dominant,  typically  flooding  even  at  ebb  tide.  The 
flooding  transported  a  large  amount  of  sediment  into  the  bay  and  navigation  channel. 
At  both  sites,  the  breach  occurred  because  the  beach  adjacent  to  the  jetty  became 
narrow  in  great  part  due  to  erosion  or  scalloping  along  the  bayside  intersection  of 
the  barrier  island  and  jetty  (Seabergh,  1999).  Kraus  and  Wamsley  (2003)  discuss 
these  and  other  locations  and  causes  of  potential  breaching  at  other  inlets  stabilized 
by  jetties  in  the  United  States. 

Most  detailed  measurements  and  numerical  models  pertain  to  breaching  of  dikes, 
which  differs  from  coastal  breaching  by  being  catastrophic  and  lasting  on  the  order 
of  hours,  and  by  the  absence  longshore  sediment  transport  that  would  tend  to  close 
the  inlet.  Visser  (1998b)  summarizes  the  literature  of  dike  breaching  in  general 
and  recent  modeling  (Steetzel  and  Visser,  1992;  Visser,  1988,  1994,  1998a)  and 
measurements  (Visser,  Vrijling,  and  Verhagen  1990;  Visser  et  al.  1995;  de  Loof  et  al. 
1996)  conducted  by  the  Dutch  research  community  (Kraak  et  al.  1994). 

Concerning  predictive  numerical  modeling,  Visser  (1994,  1998)  describes  a  five- 
stage  two-dimensional  model  of  sediment  transport  and  morphology  change  (depth 
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Fig.  3.  Breach  at  south  jetty  of  Grays  Harbor,  February  1994,  2  months  after  inception. 


and  width  of  breach)  for  dike  breaching  that  is  driven  by  the  head  difference  at  the 
dike.  The  model  was  developed  and  verified  through  laboratory  flume  measurements 
and  a  large-scale  field  campaign.  The  model  is  specific  to  dike  breaching,  which  ap¬ 
pear  to  have  different  morphologic  development  than  coastal  breaching  in  exhibiting 
an  irregular  (humped)  channel  bottom  and  a  scour  hole,  among  other  features.  How¬ 
ever,  the  stages  of  breach  development  and  associated  transport  regimes  identified 
by  Visser  (1994,  1998)  may  have  general  applicability.  Basco  and  Shin  (1999)  present 
a  one-dimensional  (elevation)  model  for  storm  breaching.  Their  model,  based  in  part 
on  the  SBEACH  dune  erosion  model  (Larson  and  Kraus,  1989;  Larson,  Kraus,  and 
Byrnes,  1990),  contains  modules  to  erode  the  dune  and  beach  across  the  barrier 
island.  Hydrodynamic  modules  are  also  included  to  furnish  wave  and  water-level 
information  across  the  calculation  domain.  Sensitivity  tests  indicated  satisfactory 
performance  that  yielded  intuitively  anticipated  results. 

From  the  above,  the  typical  coastal  breaching  process  is  as  follows.  The  ini¬ 
tiation  and  course  of  breaching  depend  on  the  cross  section  of  the  barrier  island 
or  spit,  surge  hydrograph  (or  water  level  with  respect  to  the  ocean  in  the  case  of 
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breaching  from  the  bay  or  lagoon  side),  presence  of  waves,  gross  longshore  trans¬ 
port  rate,  hydraulic  efficiency  according  to  the  location  of  the  breach  within  the 
bay  to  other  inlets,  sedimentary  composition  of  the  site,  and  presence  of  barriers 
such  as  a  jetties,  groins,  and  headlands.  In  the  absence  of  erosion-resistant  bound¬ 
aries,  coastal  breaches  appear  to  open  rapidly,  then  gradually  widen  and  deepen 
to  some  equilibrium  dimensions.  The  rapid- growth  stage  is  referred  here  as  incipi¬ 
ent  breaching,  occurring  prior  to  strong  influence  of  longshore  sediment  transport. 
Ephemeral  breaches  tend  to  open  widely  and  deepen  only  to  about  mean  sea  level, 
whereas  breaches  that  become  more  permanent  tend  to  increase  in  depth  to  that 
of  neighboring  natural  inlets  or  to  depths  of  neighboring  inlets  before  stabilization. 
If  a  breach  tends  to  become  stable,  known  empirical  relations  for  tidal  inlets  would 
then  govern  the  morphology  (see  Kraus,  2001  for  a  summary). 

Breaching  occurs  at  lower  areas  in  the  coastal  barrier,  and  the  side  slopes  of 
the  initial  breach  tend  to  be  steep,  being  eroded  by  strong  flow,  notching,  and 
collapse  in  a  type  of  avalanching.  If  erosion-resistant  materials  are  not  present,  a 
breach  typically  widens  and  deepens  rapidly  to  approach  equilibrium  dimensions 
that  depend  on  the  factors  described  above.  Models  of  coastal  breaching  are  still  at 
a  rudimentary  stage.  The  model  formulated  here  is  based  on  observations  of  general 
tendencies  of  breaches. 


3.  Analytical  Model  of  Incipient  Breaching 

A  breach  through  a  barrier  is  idealized  as  shown  in  the  rectangular  section  in  Fig.  4. 
To  derive  a  closed-form  solution,  the  governing  equations  describing  breaching  are 


Ax  Ax 


Fig.  4.  Definition  sketch  for  breaching  model. 
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developed  under  the  assumptions  that:  (a)  sediment  (sand)  volume  is  conserved; 
(b)  the  initial  condition  is  known,  (c)  sediment  transport  can  occur  at  the  bottom 
of  the  breach  and  at  its  sides;  (d)  the  breach  will  approach  equilibrium  if  external 
forces  do  not  intervene;  and  (e)  longshore  sediment  transport  is  weak  or  can  be 
neglected.  Assumptions  (b)  to  (e)  can  be  relaxed  in  extension  of  the  present  analytic 
model  or  in  a  numerical  model.  Assumption  (e),  weak  longshore  transport  relative  to 
the  breaching  process,  limits  this  model  to  describing  incipient  breaching.  Tanaka, 
Takahashi,  and  Takahashi  (1996)  and  Kraus  (1998)  present  analytical  models  of 
river  mouth  and  inlet  cross-sectional  areas  formed  under  a  balance  of  tidal  action 
that  includes  river  discharge  and  longshore  transport  under  waves.  Those  analytical 
models  that  relate  to  the  equilibrium  condition  of  an  inlet  or  river  mouth  would  be 
a  next  logical  stage  of  the  present  model  development. 

Under  the  stated  assumptions  and  by  reference  to  Fig.  4,  conservation  of  mass 


yields  the  following  equations, 

LzAx  —  QsAt 

(i) 

for  breach  width,  and, 

LxAz  =  QsAt 

(2) 

for  breach  depth,  where  L  is  length  of  the  breach  through  the  coastal  barrier,  z  and 
x  are  depth  and  width  of  the  breach,  respectively,  t  is  time,  and  Qs  and  Qb  are 
net  transport  rates  along  the  sides  (assumed  to  be  equal  on  each  side)  and  bottom 
of  the  breach,  respectively.  In  general,  the  sediment  transport  rates  Qs  and  Qb  are 
functions  of  the  acting  hydrodynamic  forces  and  breach  configuration. 

To  proceed  with  an  analytic  solution,  transport  by  surge,  or  by  flood  or  ebb  tide  is 
not  distinguished  in  forming  the  net  transport  rate  in  a  morphologic  representation, 
and  the  following  functional  relations  are  taken: 

Qs  =  Qs  _  ^  5  Qb  =  Qb  ^1  -  —  ^  (3) 

where  Qs  and  Qb  are  constant  maximum  transport  rates  that  are  not  necessarily 
equal,  and  xe  and  ze  are  values  of  the  breach  width  and  depth,  respectively,  if  the 
breach  achieved  equilibrium.  The  morphologic  parameterization  of  the  transport 
rates  by  Eq.  (3)  makes  them  time  dependent,  with  the  net  transport  rates  going 
to  zero  as  equilibrium  is  approached.  The  equilibrium  values  xe  and  ze  could  be 
specified  in  at  least  three  ways,  as  through:  (1)  experience  with  previous  breaches 
at  or  near  the  study  site,  (2)  specification  via  an  empirical  formula  such  as  given 
by  Mehta  (1976),  Graham  and  Mehta  (1981),  and  Shigemura  (1981),  or  (3)  through 
operation  of  a  hydrodynamic  model,  by  which  the  critical  shear  stress  or  similar 
diagnostic  quantity  indicating  cessation  of  transport  can  be  calculated.  The  mor¬ 
phologic  approach  as  taken  in  forming  Eq.  (3)  is  discussed  by  Kraus  (2001)  as  a 
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means  of  assuring  robust  results  though  imposition  of  morphologic  constraints  —  in 
this  situation,  of  equilibrium  values. 

Insertion  of  Eq.  (3)  into  Eqs.  (1)  and  (2)  and  taking  the  limit  yields  the  following 
coupled,  first-order  nonlinear  governing  equations  for  breaching: 


and 


x(0)  =  .x'o  >  0 


(4) 


z(0)  =  zq  >  0 


(5) 


Equations  (4)  and  (5)  cannot  be  solved  without  specifying  a  non-zero  perturbation 
of  the  barrier  (representing  an  indentation  or  “pilot  channel”  through  the  barrier), 
as  given  by  a  finite  initial  width  xq  and  initial  depth  zq.  Equation  (4)  describes  a 
one-sided  breach  such  as  constrained  by  a  jetty  in  nature  or  a  wall  in  a  flume.  If  both 
sides  of  the  breach  can  move,  then  the  value  of  Q s  in  Eq.  (4)  should  be  doubled. 


3.1.  Solution  for  xe  and  ze  — >  oo 


This  special  case  corresponds  to  very  short  elapsed  time  of  incipient  breaching,  for 
which  the  governing  Eqs.  (4)  and  (5)  reduce  to  dx/dt  =  a/z  and  dz/dt  =  b/x.  The 
solution  of  this  simplified  set  of  coupled  non-linear  equations  is, 


x  = 


a±b  Xb/a  ' 

x^a  T  (a  +  6) — - — t 


a 

CL-\-b 


(6) 


and 


z  — 


a+b  7a/b  ' 

"o  b  +  (a  +  b)  — - — t 
x0 


b 

CL-\~b 


(7) 


If  the  initial  perturbation  or  pilot  channel  is  small,  the  cross-sectional  area  of 
the  breach  is  found  to  grow  as 

/  .  Qs  +  Qb, 

xy  ~  (a  +  b)t  = - - - 1  (8) 

±j 

indicating  a  linear  increase  in  area  immediately  after  the  breach  occurs.  For  the 
special  case  of  a  =  6,  Eqs.  (6)  and  (7)  simplify  to, 

/  \  1/2 

x  =  (  x\  +  2a— tj  (9) 

and 


Zq  +  2a— t 

Xq 


z  — 


1/2 


(10) 
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which,  for  a  small  initial  perturbation,  yields  xy  ~  2 at,  a  special  case  of  Eq.  (8). 
Equations  (9)  and  (10)  make  apparent  the  significant  role  the  dimensions  of  the 
initial  perturbation  take  on  the  course  of  breach  development.  Although  Eqs.  (6)- 
(10)  govern  incipient  breaching,  they  have  limited  value  in  being  restricted  to  very 
short  elapsed  time.  More  informative  solutions  are  obtained  by  solving  Eqs.  (4)  and 
(5),  which  incorporate  modification  and  control  of  breach  growth  by  the  inclusion 
of  equilibrium. 


3.2.  Solution  with  equilibrium 


Closed-form  solution  of  Eqs.  (4)  and  (5)  is  not  found  feasible  for  arbitrary  values  of 
the  factors  a  and  6,  but,  for  the  situation  of  a  —  b  (signifying  Qs  —  Qb  —  Q)^  the 
following  solution  emerges: 

X  =  xe(l  -  f(x)e~t/T )  (11) 

and 

z  =  ze(l  -  g(z)e-^)  (12) 


in  which, 


=  Yl 
Q 


where  Ve  =  xezeL  is  the  volume  of  the  breach  at  equilibrium,  and 


/(*) 


/  1  —  ax/xe 

\1  -  OLXQ/xe 


1  /a 


(13) 


(14) 


Xp/Xe  ~  Zq/ %e 
(1  -  zp/ze)x p/xe 


fl  ZPz/Ze\1/0 
\1  -  f3z0/ze) 


(15) 

(16) 


(1  -  xo/xe)z0/ze 

Iteration  is  required  to  evaluate  Eqs.  (11)  and  (12)  because  of  the  appearance  of  the 
dependent  values  in  Eqs.  (14)  and  (16).  Equations  (4)  and  (5)  can  also  be  solved 
numerically,  as  through  a  Runge-Kutta  procedure,  and  this  was  done  to  confirm  the 
validity  of  the  closed-form  solutions  above,  as  well  as  to  perform  calculations  for 
a  ^  b. 

In  addition  to  making  basic  dependencies  of  breach  growth  explicit,  Eqs.  (11)  and 
(12)  provide  closed- form  solutions  with  which  to  check  numerical  solutions  of  the 
simultaneous,  non-linear  governing  Eqs.  (4)  and  (5).  This  simple  morphologic  model 
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indicates  that  breach  growth  is  controlled  by  seven  parameters:  xo,  zo,  xe,  ze,  Qs, 
Qb ,  and  L.  Equations  (11)— (17)  reveal  morphologic  functional  dependencies  that 
describe  the  approach  of  a  breach  toward  equilibrium.  The  quantity  r  =  xezeL/Q 
is  a  characteristic  morphologic  time  scale  governing  growth  toward  equilibrium  for 
a  given  maximum  net  transport  rate  of  sediment  removed  from  the  breach.  The 
dimensions  of  the  initial  perturbation  or  pilot  channel  of  the  barrier  island  exert 
great  control  on  the  time  development  of  the  breach  and  whether  it  will  tend  to 
widen  more  than  deepen  at  a  greater  rate,  or  vice  versa,  prior  to  approaching  equi¬ 
librium  width  and  depth.  Such  properties  of  the  solution  are  explored  in  the  next 
section. 

Equations  (11)  and  (12)  indicate  an  exponential  growth  of  a  breach  toward  equi¬ 
librium,  giving  a  more  rapid  growth  initially,  followed  by  gradual  increase  in  depth 
and  width  to  equilibrium.  The  time  behavior  of  the  morphologic  model  qualitatively 
describes  the  growth  of  breaches  observed  in  nature  and  in  the  laboratory,  whether 
induced  by  storm  surge  or  by  a  difference  in  water  level  on  the  sides  of  the  barrier 
island. 

The  volume  of  the  breach  is  V  —  xzL ,  with  x  and  z  given  by  numerical  solution 
of  Eqs.  (4)  and  (5),  respectively,  for  a  general  situation,  or  by  Eqs.  (11)  and  (12)  for 
the  special  case  Qs  —  Qb-  The  depth  of  the  breach  is  measured  from  the  top  of  the 
barrier  island  in  the  morphologic  breach  model.  For  stable  inlets,  empirical  formulas 
are  available  for  estimating  channel  cross-sectional  area  (e.g.  Jarrett  (1976)  for  large 
tidal  inlets;  Byrne,  Gammish,  and  Thomas  (1980)  for  small  tidal  inlets).  The  depth 
corresponding  to  this  channel  cross-sectional  area  is  measured  from  mean  sea  level 
to  the  bottom  of  the  breach  and  not  from  the  top  of  barrier  island.  Therefore,  in 
applications  of  predictive  expressions  for  equilibrium  breach  area  xeze)  one  must 
account  for  the  distance  from  the  top  of  the  barrier  island  to  the  elevation  of  mean 
sea  level. 


3.3.  Sensitivity  tests  of  morphologic  breach  model 

Equations  (4)  and  (5)  were  solved  numerically  for  general  cases,  after  first  confirming 
the  numerical  solution  with  the  analytical  solution  given  by  Eqs.  (11)  and  (12). 
Figures  5-7  plot  calculations  for  Qs  —  500  m3/day,  Qb  —  1,  000  m3/day,  L  —  300  m, 
and  equilibrium  width  xe  —  300  m,  and  equilibrium  depth  ze  —  5  m.  These  values  are 
considered  representative  of  the  many  small  breaches  along  the  Texas  and  Louisiana 
coasts,  as  determined  in  the  literature  review.  Note  that  if  the  crest  of  the  barrier 
island  lies,  for  example,  3  m  above  mean  sea  level,  then  the  depth  of  the  breach  below 
mean  sea  level  is  2  meters.  For  these  calculations,  the  width  of  the  pilot  channel  or 
low  section  in  the  barrier  was  specified  as  xq  =  10  m,  correspondingly,  say,  to 
walkway  or  blowout  through  the  dunes,  and  results  were  plotted  for  initial  depths 
zo  =  0.1,  1,  2,  and  3  m.  Plots  are  normalized  by  the  corresponding  equilibrium  value. 
Volume  of  the  breach  is  predicted  to  be  relatively  insensitive  to  initial  breach  depth 
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Fig.  5.  Solution  of  Eq.  4  for  breach  width,  different  values  of  initial  breach  depth. 


Fig.  6.  Solution  of  Eq.  5  for  breach  depth,  different  values  of  initial  breach  depth. 


(Fig.  7).  Growth  in  volume  toward  equilibrium  starts  faster  for  larger  initial  depth, 
but  then  the  curves  cross  at  about  30  days  elapsed  time.  The  volume  then  grows 
slightly  faster  for  smaller  initial  depth  because  of  the  larger  value  of  the  transport 
rate  at  the  bottom  than  on  the  side. 


Analytical  Model  of  Incipient  Breaching  of  Coastal  Barriers  523 


Days  after  Start  of  Breach 

Fig.  7.  Growth  of  breach  volume  for  fixed  initial  width  and  different  initial  depths  corresponding 
to  Figs.  5  and  6. 


Days  after  Start  of  Breach 

Fig.  8.  Solution  of  Eq.  (4)  for  breach  width,  different  values  of  initial  breach  width. 


The  morphologic  model  predicts  that  breach  width  grows  faster  for  smaller  ini¬ 
tial  depth.  Of  course,  breach  depth  approaches  equilibrium  faster  for  greater  initial 
depth.  Growth  of  the  dimensions  of  a  breach  depends  directly  on  the  transport  rates 
Qs  and  Qs,  and  so  general  conclusions  cannot  be  made.  However,  if  these  rates  are 
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Fig.  9.  Solution  of  Eq.  (5)  for  breach  depth,  different  values  of  initial  breach  width. 


Days  after  Start  of  Breach 


Fig.  10.  Growth  of  breach  volume  for  fixed  initial  depth  and  different  initial  widths  corresponding 
to  Figs.  8  and  9. 


of  comparable  magnitude,  as  taken  here,  then  the  trends  shown  are  qualitatively 
correct. 

Similar  calculations  are  shown  in  Figs.  8-10  for  fixed  initial  depth  zq  —  1  m  and 
initial  widths  xq  =  1,  10,  50,  and  100  m.  As  expected,  breach  width  approaches 
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equilibrium  faster  for  larger  initial  width,  whereas  breach  depth  approaches  equilib¬ 
rium  faster  for  small  initial  width.  The  volume  of  the  breach  approaches  equilibrium 
faster  for  larger  initial  width  both  because  of  the  large  initial  width  and  because 
Qb  >  Qs  h1  this  example. 


4.  Simulation  for  1980  Breach  at  Moriches  Inlet,  New  York 

Moriches  Inlet  is  a  federally  maintained  entrance  located  on  the  eastern  shore  of 
Long  Island,  New  York,  connecting  Moriches  Bay  (part  of  the  large  Great  South 
Bay)  to  the  Atlantic  Ocean  (Fig.  11).  Without  stabilization  by  structures,  inlet 
channels  on  eastern  Long  Island  meander  and  have  maximum  depth  of  about  3  m 
with  respect  to  mean  lower  low  water.  The  ocean  mean  tide  range  at  Moriches  Inlet 
is  about  1  m.  Moriches  Inlet  has  a  well-documented  history  of  opening  and  closing 
in  the  past  century  (Czerniak,  1977;  Schmeltz  et  al. ,  1982).  The  modern  inlet  was 
opened  by  a  storm  in  1953  during  construction  of  dual  jetties  that  are  spaced  245  m 
apart. 

After  inlet  opening  in  1953,  the  barrier  island  on  the  eastern  (updrift)  side  of 
Moriches  Inlet  gradually  narrowed,  primarily  due  to  erosion  at  the  jetty  and  shoreline 
intersection  on  the  back  bay.  Sorensen  and  Schmeltz  (1982)  indicate  the  adjacent 
beach  to  have  been  about  100  m  wide  in  1972  (taken  to  be  L  in  the  simulations 
described  below).  A  storm  occurring  on  January  14-16,  1980  created  a  breach  at 
the  narrowest  point  in  the  barrier  island,  about  300  m  east  of  the  east  jetty.  By 
May  1980,  the  breach  had  enlarged  to  reach  the  jetty  and  grow  to  approximately 
885  m  width  according  to  an  October  observation  (see  Fig.  12  for  a  September 


Fig.  11.  Location  map  for  Moriches  Inlet,  Long  Island,  New  York. 


526  N.  C.  Kraus 


Fig.  12.  Breach  at  Moriches  Inlet,  New  York,  September  21,  1980.  Photograph  taken  9  months  after 
the  January  1980  breach.  Note  east  jetty  became  isolated  from  land  because  of  the  breach. 


Table  1.  Dimensions  of  1980  breach  at  Moriches  Inlet,  New  York  (data  from 
Schmeltz  et  a/.,  1982,  and  Sorenson  and  Schmeltz,  1982). 


Date,  1980 

Estimated 
Elapsed  Time, 
Days 

Breach  Width, 
m 

Breach  Depth  from 
Assumed  Elevation  of 
Barrier  Island  Crest,  m 

15  January 

0 

30  (assumed  here  for  :co) 

1  (assumed  here  for  zq) 

16  January 

1 

90 

3.6 

20  January 

5 

215 

4 

“May” 

90 

885 

7 

observation),  more  than  three  times  the  width  between  the  two  jetties,  and  with  a 
depth  of  3  m  relative  to  mean  low  water  (Sorensen  and  Schmeltz,  1982). 

Table  1  compiles  available  information  about  the  breach  dimensions.  For  the 
tabulation,  it  was  assumed  that  the  mean  elevation  of  the  barrier  island  was  3  m 
above  mean  low  water.  This  value  was  added  to  the  depths  of  the  breach  reported 
by  Sorensen  and  Schmeltz  (1982).  In  addition,  an  initial  width  of  30  m  and  initial 
depth  (local  lowering  in  the  barrier  island)  of  1  m  were  assumed  for  application  of 
the  breach  morphology  model.  Based  on  additional  photographic  evidence,  it  was 
concluded  for  this  analysis  that  near-maximum  width  had  been  achieved  sometime 
before  the  May  observation.  This  value  was  set  in  the  model  at  90  days  elapsed 
time  after  opening  of  the  breach.  Construction  began  in  October  1980  to  close  the 
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breach,  for  which  900,000  m3  were  dredged  from  the  bay  and  trucked  from  a  quarry 
in  equal  amounts.  This  volume  provides  a  check  for  the  model. 

The  breach  morphology  model  was  established  for  the  area  with  the  initial  di¬ 
mensions  as  in  Table  1.  Trial  and  error  readily  gave  a  visual  good  fit  to  the  data 
with  base  rates  of  Qs  —  Qb  —  30,000  m3/day  for  the  maximum  transport  rates  at 
the  inception  of  breaching.  Such  rates  are  well  within  the  daily  production  capacity 
of,  for  example,  0.5-  and  0.6-m  (20-  and  24-inch)  diameter  dredges  operating  in  shal¬ 
low  water  (Turner,  1996)  and,  therefore,  are  considered  readily  attainable  by  storm 
surge  combined  with  tidal  exchange.  To  explore  sensitivity  of  the  model  to  variations 
in  maximum  transport  rates,  the  base  rates  were  multiplied  by  2  and  0.5  to  define 
plausible  upper  and  lower  bounds,  respectively.  Figures  13-15  plot  calculations  with 
the  observed  width,  depth,  and  estimated  volume  of  the  breach,  and  the  results  with 
upper  and  lower  limits  of  assumed  maximum  transport  are  depicted  with  dashed 
lines.  The  volume  at  different  times  was  estimated  from  the  dimensions  as  given 
in  Table  1  and  the  reported  fill  volume,  assumed  to  occupy  the  same  approximate 
equilibrium  volume  in  May. 

The  breach  morphology  model  captures  the  qualitative  behavior  of  the  breach 
evolution,  with  rapid  initial  growth  followed  by  gradual  approach  to  equilibrium. 
Width,  depth,  and  volume  of  the  breach  are  in  quantitative  agreement  with  ob¬ 
servations  for  all  available  data  points  through  time.  Calculations  with  the  lower 
estimate  of  the  maximum  transport  rate  deviate  more  from  the  observations  than 


Fig.  13.  Simulation  of  growth  of  breach  width,  Moriches  Inlet. 
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Fig.  14.  Simulation  of  growth  of  breach  depth,  Moriches  Inlet. 
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Fig.  15.  Simulation  of  growth  of  breach  volume,  Moriches  Inlet. 
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those  for  the  base  rate  and  double  the  base  rate.  The  characteristic  morphology 
growth  time  is  r  =  Ve/Q  =  900,  000/30,  000  =  30  days. 


5.  Discussion  and  Conclusions 

This  paper  has  presented  an  analytical  morphologic  model  of  incipient  breaching 
of  alluvial  barriers.  Closed-form  solutions  were  obtained  through  simplifying  as¬ 
sumptions  about  breach  configuration  and  approach  to  equilibrium.  Processes  as¬ 
sociated  with  the  driving  hydrodynamics  of  sediment  transport  and  bank  erosion 
and  collapse  are  represented  implicitly  through  specification  of  the  maximum  sed¬ 
iment  transport  rates  at  the  bed  and  the  sides  of  the  breach  at  initial  breaching, 
and  by  the  equilibrium  values  of  the  breach  width  and  depth.  Additional  or  micro¬ 
scale  physical  processes  could  be  represented  in  numerical  solution  of  the  governing 
Eqs.  (4)  and  (5),  for  which  time-dependent  sediment  transport  rates  and  bank  fail¬ 
ure  could  be  calculated  while  maintaining  a  known,  geometry.  Seven  variables  are 
found  to  control  breach  growth  in  such  a  macroscale  description:  initial  width  and 
depth  of  the  breach,  equilibrium  width  and  depth  of  the  breach,  width  of  the  bar¬ 
rier  island,  and  maximum  or  initial  net  sediment  transport  rates  at  the  bottom  and 
sides  of  the  breach.  Breach  growth  follows  an  exponential  behavior  governed  by  a 
characteristic  time  scale  r  =  Ve/Q,  where  Ve  is  the  volume  of  the  breach  at  equilib¬ 
rium,  and  Q  is  a  representative  maximum  net  sediment  transport  rate  through  the 
breach. 

Sensitivity  testing  of  the  morphologic  model  and  comparison  with  observations  of 
the  1980  breach  at  Moriches  Inlet,  New  York,  demonstrated  validity  of  the  model  in 
capturing  the  general  qualitative  and  quantitative  features  of  coastal  barrier  breach¬ 
ing.  Data-based  refinements  can  be  made  by  the  inclusion  of  empirical  predictive 
formulas  for  the  equilibrium  depth  and  width.  Most  of  these  extensions  would  re¬ 
quire  information  on  the  hydrodynamics  and  sediment  transport  at  the  breach, 
which  would  need  to  be  supplied  by  hydrodynamic  and  sediment  transport  models 
coupled  to  the  morphologic  breaching  model.  However,  estimates  of  the  basic  gov¬ 
erning  variables  drawn  from  engineering  judgment  and  observations  at  the  site  and 
at  nearby  inlets  may  be  suitable  for  obtaining  a  reconnaissance-level  prediction  of 
coastal  breach  growth  without  recourse  to  sophisticated  models. 
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